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Abstract 

The multiplicity distribution and the two-particle Bose-Einstein correlations at fixed multiplic- 
ities observed in pp collisions at y/s = 7 TeV by the ALICE Collaboration are analyzed by the 
formulae obtained in the quantum optical approach. The chaoticity parameters in the inclusive and 
semi- inclusive events are estimated from the analysis. Multiplicity or fcr dependence of longitudinal 
and transverse source radii are also estimated. 
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I. INTRODUCTION 



In high energy nucleus-nucleus or hadron-hadron collisions, Bose-Einstein correlations of 
identical particles are considered as one of the possible measures for the space-time domain 
where those particles are produced. As the colliding energy increases up to the TeV re- 
gion, Bose-Einstein correlations of identical particles can be precisely analyzed in the semi- 
inclusive events or at fixed multiplicities to estimate the production domain of identical 
particles. 

Up to the present, most of theoretical approaches to identical particle correlations in the 
semi-inclusive events are investigated in the case of purely chaotic field [l-4]. 

One of the theoretical approaches to the Bose-Einstein correlations is made on the anal- 



ogy of the quantum optics 



, |6[, where two types of sources, chaotic and coherent sources 

achs formula [^j which 
to find higher order 



are introduced. A diagrammatical method, based on the Glauber- 
is derived from the single-mode laser optics, has been proposed 3] 
Bose-Einstein correlation (BEC) functions in the inclusive events. In Ref . 8] , the generating 
functional (GF) for the momentum densities is formulated, and a diagrammatic representa- 
tion for the cumulants is proposed. Identical particle correlations in the semi-inclusive events 

n n 

are formulated in [9[. It is shown [lOj that the multiplicity distribution (MD) in the quantum 
optical (QO) approach is approximately given by the Glauber-Lachs formula jol. Illj]. 

In general, in the formulation of two-particle BEC function at fixed multiplicities or in the 
semi-inclusive events, two-particle and one-particle momentum densities in the semi-inclusive 
events and the MD are contained. If the two-particle momentum density is integrated over 
one of the momenta, the one-particle momentum density is obtained. If the one-particle 
density is integrated over, the MD is obtained 12j . Therefore, even in the observed MD, 



information on the source radii would be included. 

In the present paper, the MD and the two-particle Bose-Einstein correlations in the semi- 
inclusive events observed by the ALICE Collaboration are analyzed by the formulae obtained 
in the QO approach. The chaoticity parameter p- m in the inclusive events is estimated from 
the analysis of the observed charged MD. The chaoticity parameter p sm in the semi- inclusive 
events is estimated from the analysis of two-particle Bose-Einstein correlations. Longitudinal 
and transverse source radii at fixed multiplicities are also estimated from the data. 



II. MOMENTUM DENSITIES IN THE SEMI-INCLUSIVE EVENTS 



In the QO approach, the n-particle momentum density in the semi-inclusive events is 
defined by, 



p n (Pi,---,Pn) = co(\f(Pi)\ 2 ■ ■ ■ \ f(Pn)\ 2 ) a , n = 1,2,... 

M 

f(P) = ^°i&(p) + fc(p), 



(1) 



i=l 



where Co is a normalization factor. In Eq.flT]), 4>i(p) and f c {p) are amplitudes of the ith 
chaotic source and a coherent source, respectively, and a, is a random complex number 
attached to the ith chaotic source. In addition, M is the number of independent chaotic 
sources^], which is regarded to be infinite in the present paper. 

In Eq.flT]), parenthesis (F) a denotes an average of F over the random number a« with a 



Gaussian weight 



6] 



M I |2 

W-n(^//exp[-^K<,)K 



(2) 



One-particle and two-particle momentum densities are respectively given as, 

Pi(Pi) = c (\f(pi)\ 2 } a = c [r(pi,pi) + c(pi,pi)], (3) 

P2{P1,P2) = C (\f(pi)f(p2)\ 2 }a 

= c o{p(Pi)p(P2) + \r(pi,p 2 )\ 2 + 2Re[r(p 1 ,p 2 )c(p 2 ,p 1 ))}, (4) 

where r(pi,p 2 ) is a correlation caused by the chaotic sources and c(pi,p2) is a correlation 
by the coherent source, 

M 

r(pi,P2) = ^Ai^i(pi)0j(p2), c(px,p 2 ) = /c(pi)/*(P2). (5) 
1=1 

The generating functional (GF) of momentum densities in the semi-inclusive events is 
defined by the following equation [l2|, 

d 3 p 1 d 3 p n 



n=l ' J J 



p n (pi, . . .,p n )h(pi) ■ ■■h{p n )- 



E„ 



(6) 



From Eqs.fllJ, (J2J) and flBJ), the GF is written as, 



Z sm [h(p)] = c (exp 




fip) I 2 ^(P) 



(7) 



where an additional constant Z sm [h(p) = 0](= Co) is added to the right hand side of Eq.Q. 

Inversely, the n-particle momentum density in the semi-inclusive events is given from the 
GF as 

5 n Z sm [h(p)} 



p n (Pl, ...,p n ) = E 1 ---E ri 



The nth order cumulant is given by 



g n (pi, . ..,p n ) = E x ■ ■ -E n 



8h(pi) ■ ■■5h(p n ) 



S n \nZ sm [h(p)} 



(8) 



h(p)=0 



8h(pi) ■ --Shipr, 



(9) 



h(p)=0 
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From Eqs.flEl) and (|9]), we have an iteration relation for the momentum densities [9], 
Pi(Pi) = co^i(pi), 

p n (p U . . . ,p n ) = gx(pi)p n -l(P2, ...,p n ) + C g n (Pl, ■ ■ -,Pn) 

n-2 

, Pj 1 , . . . , Pj i 



(10) 



i=l 



The second summation on the right hand side of Eq. ( flOl indicates that all possible combi- 
nations of (ji, . . . ,ji) and (ji+i, . . . ,j n -i) are taken from (2, 3, ... , n). 

In order to calculate momentum densities at fixed multiplicities, following equations are 
defined for n = 1, 2, • • • and k — 0, 1, • • • , n — 1; 

d 3 p k+1 d 3 p n 



Pn\Pu---,Pk) : 
g { n k \pi,...,p k ) 
Then the MD is given by 



1 

(n-k)\ 
1 

(n - k)\ 



p n (pi, . . . ,p k ,Pk+l, ■ ■ ■ ,p n y 



E, 



fc+1 



Er, 



,d 3 p k+1 d 3 p n 

g n {Pl,---,Pk,Pk+l,---,Pn)-F l F — ' I 11 ) 



E k +i 



P(0) = c , 



P(n) = p. 



(0) 



(n-k)\ 



P ( n ] (Pl,---,Pk) 



d 3 pi d 3 p k 
E\ E k 



(12) 



From Eqs.flHID, (ffU) and we have 

1 n 

P(n) = -J2j9? ] P(n-j) 



(13) 



n j-1 

P ( u\pi) = Y.^f\pi)P{n-j), gf\ Pl ) = Rj{PiiPi) +^ T y-w(^^i)' ( 14 ) 

3=1 1=0 
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n-l 



Pn\PuP2) = ^2( n -f)9j (Pl)P n -j(P2) + (Pi, Pi) P(n ~ j) , 

J=l J=2 

9j 2) (Pl,P2) = ^2Rj(pi,p 2 )Rj-l(p2,Pl) 

1=1 

+ ^^{T m ^ m (p 1 ,p 2 )Rj-i-i(p2,Pi) + Rj~i-i(puP2)T m ,i- m (P2,Pi)}, (15) 

Z=0 m=0 

where, with Ro(ki, k 2 ) = Ui5 3 (ki — fc 2 ), 

<i 3 fc . (<n f f f „ . , d 3 A; 



Af » = /// Rj(k, kf-±, Af , = 111 TJ^ft (16) 
Rj(pi,P2) = J J J r(pi, k)Rj-i{k,p 2 )-^-, 

Tj,i(pi,P2)= [ ■ ■ ■ [ R j (pi,k 1 )c(k 1 ,k 2 )Ri(k 2 ,P2)^-^^-^ L . (17) 
J J Wi u 2 

IfEq.(H5D for pf( Pl 

,p 2 ) is integrated over p 2 , Eq.dHD for pn \pi) is obtained. In addition, 
if Eq. CHj) is integrated, Eq. (IT5|) . the recurrence equation for the MD, is obtained. 



III. PARAMETRIZATION 



In the followings, variables are changed from (ph,Pit) to (yi,Pir) with rapidity y\ 
defined by y± = tanh~ 1 (piz / /£'i). Correlations r(pi,p 2 ) and c(pi,p 2 ) are both assumed to be 
real, and are parametrized as [9], 

r{yi,PiT, V2, P2t) = P S mV p(yi, Pvr)p(V2, P2t) /(Ay, Ap 1T ), 

c(yi,PlT] 2/2, P2T) = (1 - Psm) V P(VU PlT)p(y2, Ptt) , (18) 

p(Vi,Pit) = K)-^^exp[-a?/i 2 -/3pi T 2 ], 
I(Ay, Ap T ) = exp[- lL {Ay) 2 - 7t(A Pt ) 2 ], (19) 

where Ay = y 2 — yi and Apx = P2T — Pit- Six parameters, p sm , (no), a, (3, 7^ and 
7t are included in Eqs.( lT8"j) and ( jTSJ) . All of them are assumed to be positive constant 
at present: Parameter p sm denotes the chaoticity parameter in the semi- inclusive events 
with < p sm < 1. Parameter (no) is related to an average multiplicity. In the limit of 
Psm = 0, the MD defined by Eq. (|T2|) or Eq. (|T3|) becomes the Poisson distribution with an 
average (n ) . Parameters a and (3 are related to the width of rapidity and px distributions, 
respectively. The longitudinal momentum transfer squared, qf = (E 1 — E 2 ) 2 — (pu ~P2l) 2 , 
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is approximately written as qf on „ ~ 2(mT 2 )(cosh Ay — 1) with the average of transverse mass 
squared (my 2 )- Then, g 2 ong w {mr 2 ){Ay) 2 for \Ay\ << 1. Therefore, (mx 2 ) is roughly 

equal to the longitudinal source radius. Similarly, ^/yr denotes the transverse source radius. 

Then, the jth order correlation, Rj(pi,p 2 ), of chaotic component is written in the follow- 
ing form, 



Rj(yi,PiT,y2,Pzr) = Njexpl-Ajtyt 2 + y 2 ) + 2Cfyxy^ 

x exp[-^(p 1T 2 + p 2T 2 ) + 2W jPlT p 2T }, (20) 



where 



W x = 7 T , W j+1 



AT i \ V®P AT p sm (n Q )x/af3 

Nl = Psm(^o) T77T-, = ; -/Vi- 

^sm\ 0/^, ,+1 ^ +Al (^. + f/l ) 



Recurrence equations for Aj, Cj, Uj, Wj and Nj can be solved |3|,HQ( , and 

r 2 - n 1 + (ri/r 2 ) J _ _ (r 1 /r 2 ) j/2 

J ~ 2 l-(n/r a )J' J_lr2 riJ l-(r!/r 2 )i 

„ _ * 2 - «i i + (h/t 2 y w (h/t 2 y/ 2 



2 l-(ti/f a y J w "l-faAa)'' 
= y^t 2 / Psm (noW /2 P y f 1 - (n/r 2 ) i V2 1 - (tx/t 2 ) 
J 7r 3 /2 \ r2 i/2 t2 ) \l-( ri /r 2 y) l-ih/t^' { ' 



where 



a 



+ 2y L - \/a 2 + 4<xy L _ a + 2y L + \/ a 2 + Aa^ L 



n = i , r 2 



2 7 2 

P + 2 7t - x/fi 2 + A(3 lT f3 + 2 lT + x/fi 2 + Af3 l7 



(22) 
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IV. ANALYTICAL FORMULA FOR MULTIPLICITY DISTRIBUTION 



From Eq. (|2ip . the recurrence equation for the MD, Eq. (|13p . is written as 



1 n 

P(n) = [A^ + j^j-i] P(n - j), 
U i=i 



n = l,2,..., 



A 



(R) 



in- 2 



(23) 



where 



A = a/1 - ri/r 2 (l -t 1 /t 2 )(l - Pea) (no). 



(24) 



As can be seen from the above equations, the MD contains four parameters, p sm , (no), 
hh = 1l/oc and h T = j T /P- 

The analytical formula for the MD can be obtained approximately from Eg. (1231) [10]. The 
generating function for P(n) is defined by 

oo 

II(z) = ^P{n)z n . 

3=0 

Then, the differential equation for II (2) is obtained from Eq.( l23l) as, 

d 

ll{z) = 

3=± 

It is shown from Eq. (1221) that < r x /r 2 < 1 and < 1 1 /t 2 < 1- Therefore, for ri/r 2 , t\jti << 
1 or j » 1, Eqs. fl23|) can be approximated by, 



, oo 

A nW = ^(Af) +j Af_>)^n(,). 



(25) 



Then, Eq.(l25l) can be rewritten as, 

dU 



5 + M^P 1 



.1 — £2 ' U ^ dz L 1 — £z. 
With the boundary condition 11(1) = 1, we obtain 

A) " 1) 



1-^— -(z-l)| exp 



:i-o 2 i-^(.-i) 



The average multiplicity (n) is given from H(z) as 



(n> 



+ 



An 



(26) 



2 =i 1-^ ' (i-0 2 ' 

From Eq. (|24|) . parameter £ is proportional top sm (nn), whereas A is to (1— p sm )(no). There- 
fore, from Eqs. fllSp and ( I19p . we can identify the first term on the right hand side of Eq. ()26p 
as a contribution of the chaotic sources and the second term that of the coherent source to 
the average multiplicity (n) ; 



Pm)(n), 



1_£ — (l_ e) 2 v 

where p; n denotes the chaoticity parameter in the inclusive events. 
Then, the generation function is given as 

(z-1) 



(27) 



U(z) = {l -pi n (n)(z - 1)} exp (l-p in )(n) 



1 -Pin{n)(z- 1) 



(28) 



The MD is given from the GF as 
1 d n 



n n ) = -r-r- n M 

v ' n\dz n v ' z=o 
(Pm(n)) n 
" (1+Bn(n))" +1 



exp 



(1 -p in )(n)] /(l -p in )(n) 



1 + Pin(n) 



n = 0,1, 



(29) 



where L„{x) denotes the Laguerre polynomial. Equation f[2"9l is called the Glauber-Lachs 



formula 15 
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111 ]. The KNO scaling function of the Glauber-Lachs formula is given by, 



1 

— exp 

Pin 



Z + l- Pn 
Pin 



tat— Vi 1 "Pin)*), 



(30) 



where Iq(z) is the modified Bessel function. 



V. TWO-PARTICLE BEC FUNCTION 

The two-particle Bose-Einstein correlation (BEC) function Cn (Ay, Ap T ) at fixed mul- 
tiplicity n is defined as 

C {2 \Ay Ap T ) = JJJpn\yi,PiT,yi + Ay, Pit + Ap T )d yi d 2 p lT 

n ~ 1 III Pn\yuPiT)pn\yi + Ay, Pit + Ap T )dyid 2 pi T 

Each term in Pn(pi,P2) or p„ (pi)pn (pi) in Eq. (l3Tj) takes the following form, 

F{yi,Pvr,y2,Pzr) = C exp[-A uyi 2 - A 22 y 2 2 + 2A 12 y 1 y 2 ] 

x exp[-U u piT 2 - U 22 p 2T 2 + 2Ui 2 piTPzr], (32) 



where 2/2 = 2/1 + Ay, Vrr = Pit + Ap T . C , An, A u , A 2 2, U u , Uu and U22 are some 
constants. 

In the observed Bose-Einstein correlations by the ALICE Collaboration 13( , the average 
kx of the two particle transverse momenta is restricted in a small region. In order to take 
the constraint into account, we parametrize, p\x = k T — Apx/2, and Pit = &t + 
Using the approximation, fey • Apx = {Pit 2 — P2T 2 )/2 ~ 0, we rewrite the terms on the 
transverse momenta in Eq. fl32|) as follows 



-UuPxt 2 - U22P2T 2 + 2U 12 PlTP2T 

= -(U n + U22 - 2U l2 )k T 2 - {Uu + U 22 + 2U l2 )A PT 2 /A. 

After the integral for p 1T in Eq.( ET|) is replaced by that for k T , the condition that \k T \ (=k T ) 
is approximately constant, \kr\ ~ c, is introduced into Eq.(|3TT); 



G(Ay, Ap T ) 



y/A n + A 22 - 2A 



F(yi,PiT,yi + Ay,p lT + Ap T )dy l 

Ay 2 



fc T |«C 



: exp 



12 



A11A22 — A\2 2 



An + A 



22 



2A 



12 



x exp 



-(U u + U22 + 2U 12 )Ap T 2 /A - (U u + U22 - 2U 12 )k T ' 



(33) 



The width of \k T \ is neglected in Eq. fl33l . because it is canceled out in Eq.flBT 



VI. ANALYSIS OF EXPERIMENTAL DATA 

The ALICE collaboration reported the observed charged multiplicity distribution [l^] 
within the pseudo-rapidity range |?7| < 1.0 and the observed two-particle Bose-Einstein 
correlations [13] within \r]\ < 1.2 in pp collisions at -y/i = 7 TeV. We analyze the data 
on Bose-Einstein correlations within \r]\ < 1.2, using the parameters estimated from the 
observed MD within \rj\ < 1.0. 

In order to analyze the observed MD by Eq. (l23|) . we extract data of even prongs and 
renormalized those data to 1, because Eq. (123]) is for identical particles. The result on 
the minimum chi-squared fit is Xmm/ n °f = 81.818/(35 — 4) at p sm = 9.15, (n ) = 0.863, 
hi = 9.36 x 10 -5 and Iit = 9.34 x 10~ 5 . Estimated values of 7^ and 7t become so small that 
we cannot fit observed Bose-Einstein correlations for longitudinal and sideward directions. 

Then, we use Eq. (!28|) . the KNO-scaling function for the Glauber-Lachs formula, because 
we can use the published data without modification, and identify the scaling variable z c ^ = 

9 



TABLE I. Estimated parameters in the analysis of charged MD observed in pp collisions at yfs = 7 
TeV 



(rich) Pin Xmin/ nof 

12.11 ±0.09 0.654 ±0.021 235.3/(65-2) 



rich/ (n c h) for charged particles with z = n/(n) for identical particles (for example negatively 
charged particles), where (n ch ) = 2{n). 

We analyze the MD observed in pp collisions at <Js = 7 TeV by P(n) = 4>(z) / (n c h) , where 
4>(z) is given by Eq. fl30l) . The estimated parameters in the analysis are shown in Table [B 
The observed MD is compared with P{n) calculated by Eq. (!30|) in FigJTJ 

1.0x1 (T 1 r i 



C 

v x 

3 

-e- 

A 

0- 1.0x1(T 2 - 



1.0x10" 3 



AUCE pp7TeV 

hl<i 



X 



Data 
■GL 




1.0x10- 

0.0 1.0 2.0 3.0 4.0 5.0 6.0 

z=n/<n> 



FIG. 1. Analysis of charged MD in pp collisions at a/s = 7 TeV by Eq. ([30]) . 

As can be seen from Eqs.( 124"j) and ( E7I) . the estimated values of (n c h) (or (n) = (n ch )/2) 
and pi n gives two constraints on the parameters, p sm , (no), hi and hr- When we analyze the 
observed data on two-particle Bose-Einstein correlations with 4 parameters, a, fl, hi and 
Ht, values of Hl = and Ht = Jt/P sometimes show unstable behaviors. Therefore we 
fix at a = 0.25 and /3 = 50. Then, we analyze the data with two parameters, and hr. 

The parameters p sm and (no) are given by the following equations, 

U2 ( v _ Aof l+ U2 

Psm ~ Aoh+Ui {n ° } ~ hh ■ 
In the above equations, ^ and A Q depend only on p [n and (n), whereas fi and /2 depend 
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only on 7l/« and Jt/P] 

£=Phx(n)/(l+Pin(n)), A = (l-p in )(n)/(l+ Pin (n)) 2 , 

fi = (1 - VnF2)(i - VtJh) 2 , h = ^l-ri/r 2 (l - ti/t 2 ). 

The two-particle Bose-Einstein correlations in collisions are measured in longitudinal, 
sideward and outward directions. In pp collisions, data would be almost the same in g si d e and 
g out directions. Therefore, the data in gi ong and q^e directions are analyzed. We compare 
the calculated results for n — 4, 10 and 17 with the data for n c h from 1 to 11, from 17 to 22 
and from 29 to 34, respectively. 

In the analysis of two-particle Bose-Einstein correlations, we use the identity Ap? 2 = 
Qside 2 + Qout 2 and the following approximate relations, 

ftong = 2(m T 2 )(cosh Ay - 1), 

Ay = ln(a±Vrf=T), a =-^4 + 1, (34) 

2{rriT ) 

where (m T 2 ) = 0.14 2 + 0.2 2 = 0.0596 (GeV/c) 2 . 

(2) 

Data on C n ((Ziong) are taken with the fcr range (/jti,&T2) GeV/c, |g S ide| < 0.16 GeV/c 
and |<7out| < 0.16GeV/c. In order to take these conditions into account, we use fc^ 2 = 
(k T i 2 + k T2 2 )/2, Ap T 2 = g side 2 + g out 2 = 2 + 0.16 2 = 0.0256 (GeV/c) 2 in Eq.(T55]). Data 
on Cn (g S ide) are taken with the k T range, (^1,^2) GeV/c, l^iongl < 0.16 GeV/c and 
|?out| < 0.16GeV/c. In this case, we use kx 2 = (kxi 2 + A>T2 2 )/2 and gi ong 2 = Qout 2 = 0.0128 
(GeV/c) 2 in Eqs.flSD and (|3l|) . 

The results on Cn^ (ftong) are shown in Table [III and FigfJl As the data are taken with kx 
range (0.2,0.3) GeV/c, |g sid e| < 0.16 GeV/c and |g out | < 0.16 GeV/c, we put k T 2 = 0.065 
(GeV/c) 2 and Apr 2 = 0.0256 (GeV/c) 2 in Eq.(l33l. In the minimum chi-squared fit, data 
points with g long < 0.6 GeV/c (15 points) are used. All the values of parameters, p sra , (no), 
hi and h T estimated in the analysis of Cn \q\ ng) increase as multiplicity n increases. It 
should be noted that the MD defined by Eq. (123]) is contained in the formula of the two- 
particle BEC function given by Eq. (13 1 1) . and it calculated from Eq. (l23l with parameters 
shown in Table Hi] becomes broader as multiplicity n increases. 

The results on Cn (<7side) are shown in Table IIHI and FigJSJ As the data are taken with 
k T range (0.2,0.3) GeV/c, \q long \ < 0.16 GeV/c and |g side | < 0.16GeV/c, we put k T 2 = 0.065 
and Ap T 2 = 0.0256 in Eq.f l33]) . In the minimum chi-squared fit, data points with q long < 1 
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TABLE II. Estimated parameters in C\ (oWf 


1 — 

r) observed in pp collisions [13 


at yjl = 7 TeV. 


parameters 


n = 4 


n = 10 


n = 17 




0.969 ± 0.004 


0.975 ± 0.003 


0.978 ± 0.003 




3.78 ±0.10 


4.93 ±0.12 


6.07 ±0.12 


h L 


5.13 ±0.29 


6.56 ±0.34 


8.29 ±0.38 


hr 


0.351 ±0.014 


0.543 ±0.013 


0.698 ± 0.001 


Xmin/ nof 


99.1/(15 - 2) 


84.0/(15 - 2) 


34.8/(15 - 2) 



3.0 



2.5 



2.0 



1.5 



1.0 



ALICE pp7TeV 

k T (0.2,0.3) c nch(1-11)+1 

BEC(n=4)+1 

nch(17-22)+0.5 

X 4 '°™ a nch(29-34) 
BEC(n=17) 



^ ~ » a a o 

ffl $ 

<>i,\j + « * 

>--<.-<>-*- 6 - -*- 



a 

0.5 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

q, (GeV/c) 

long 



(2) 

FIG. 2. Analysis of C% '{qi orig ) 

in collisions at yfs = 7 TeV by Eq.([3T 



GeV/c (25 points) are used. All the values of parameters estimated in the analysis increase, 
as multiplicity n increases. The MD calculated from Eq. (j23l) with parameters shown in 
Table II I II also becomes broader as multiplicity n increases. 

The results on the kx dependence of (qWg) are shown in Table HVl and FigJU All the 
values of parameters increase as kx increases. In this case, the MD calculated from Eq. fl23|) 
with parameters shown in Table IIVI becomes narrower as kx increases. 

The results on fey-dependence of the Bose-Einstein correlations in the sideward direction 
are shown in Table |V] and in Figj5j All the values of parameters decrease as kx increases. 
The MD calculated from Eq. f l23|) with parameters shown in Table IIVI becomes narrower as 

12 



TABLE III. 


(2) 

Estimated parameters in C n 


(fede) observed in pp collisions at •> 


1 1 

/i = 7 TeV [LJ. 


parameters 


n = 4 


n = 10 


n = 17 




0.963 ± 0.004 


0.971 ±0.003 


0.976 ± 0.003 


(no) 


3.14 ±0.05 


4.53 ±0.10 


5.79 ±0.1 


h L 


6.76 ±0.19 


10.8 ±0.3 


14.4 ±0.5 




0.106 ±0.007 


0.241 ±0.016 


0.357 ±0.022 


Xmin/ nof 


142.1/(25 - 2) 


84.0/(25 - 2) 


34.8/(25 - 2) 



3.0 



2.5 



2.0 



1.5 



1.0 



ALICE pp 7 TeV 

k (0.2,0.3) ° " ch ( 1 -") +1 
T BEC(n=4)+1 



C 4 (2) (q )+1 
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FIG. 3. Analysis of C n (q s ide) i n PP collisions at yfs = 7 TeV by Eq. ([3T 



kr increases. 

The source radii, Ri for the longitudinal direction and R s for the sideward direction are 
defined by the following equations, 



Ri 



m T 2 ), R s 



It- 



The multiplicity n dependences of source radii calculated from jl and 7t, which are 
estimated in the analysis of Cn\q\ ons ), are shown in Figj6j Source radii Ri and R s increase 
as multiplicity n increases. Values of Ri and R s at each n coincide within the error bars. 

The n dependences of source radii calculated from ^ L and 7^ in the analysis of Cn\q S ide) 
are shown in FigJTJ Both radii increase with multiplicity n. At each n, Ri is greater than 
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FIG. 4. Analysis of C} (g s ide) i n PP collisions at -y/s = 7 TeV by Eq. (f3Tj) , 



TABLE V. 


(2) 

Estimated parameters in C\ (q s idc) 


observed in pp collisions at 


v/i = 7 TeV |l4]. 


ClO^side) 


k T = (02,03)GeV 


k T = (04,05)GeV 


k T = (06,07)GeV 


Psm 


0.971 ± 0.003 


0.967 ± 0.004 


0.959 ± 0.005 


(no) 


4.35 ±0.11 


3.66 ±0.11 


2.67 ±0.12 


h L 


10.3 ±0.34 


8.45 ± 047 


4.75 ±0.52 
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0.223 ± 0.018 


0.151 ±0.014 


0.0894 ± 0.0089 
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The &>r dependences of source radii Ri and i? s calculated from ^ L and 7t in the analysis 
of Cjg^giong) are shown in FigJHl Both radii decrease as k T increases. The source radius Ri 
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FIG. 5. Analysis of Cn ((?side) m PP collisions at y/s = 7 TeV by Eq. ([3T 
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FIG. 6. Source radii Ri and R s in the analysis of C\ wiong), which are calculated from 7^ and 
7t shown in Table HU 



is equal to R s within the error bars at each kx- 

The kr dependences of source radii Ri and R s calculated from 7^ and 7-7- in the analysis 
of Cio^(fZside) are shown in Figj9j Both radii decrease as k? increases. Ri is greater than R s 
at each hp. 



VII. SUMMARY AND DISCUSSIONS 

The charged multiplicity distribution in \q\ < 1 and the two-particle Bose-Einstein cor- 
relations at fixed multiplicities in \r]\ < 1.2 observed in pp collisions at s/s = 7 TeV by the 
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FIG. 7. Source radii Ri and R s in the analysis of CA ((/side), which are calculated from 7x and 
7t shown in Table IHIi 
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FIG. 8. Source radii Ri and R s in the analysis of C 10 (q\ ng)- The radii are calculated from 7^ 
and jt shown in Table HVl 



ALICE Collaboration are analyzed by the formulae obtained in the QO approach. At first, 
the observed charged MD is analyzed by the KNO scaling function of the Glauber-Lach 
formula. It is approximately derived from the recurrence equation for the MD, and con- 
tains two parameters, the chaoticity parameter p in in the inclusive events and the average 
multiplicity (n) of negatively charged particles. Those two parameters give two constraints 
on the chaoticity parameter p sm in the semi-inclusive events, the normalization factor (no), 
hi = ji/a and hx = 7t/ ' fl- Parameters a and (3 are related to the width of rapidity and 
Pt distributions, respectively. Those are fixed at a = 0.25 and (3 = 50. The parameter 7x 
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FIG. 9. Source radii Ri and R s in the analysis of C} (g s ide)- The radii are calculated from 72, 
and 7t shown in Table PVl 



is related to the longitudinal source radius and 7^ is related to the transverse radius. Using 
parameters hi (or 7^) and hr (or 77-), we can analyze the data on two-particle Bose-Einstein 
correlations for the longitudinal direction and the sideward (or outward) direction. As the 
dada are taken in pp collisions, we neglect the difference between the correlations in the 
sideward direction and those in the outward direction. 

In our formulation of the BEC function, restrictions on kx, 9out> <?side and so forth are 
taken into account. Therefore, even in the analysis of two-particle Bose-Einstein correlations 
in the longitudinal direction, we can estimated the source radius R s in addition to the source 
radius Ri. 

From the analysis of the Bose-Einstein correlations in the longitudinal direction, esti- 
mated values of Ri and R s coincide with each other at each multiplicity n or kx- On the 
other hand, from the analysis of the Bose-Einstein correlations in the sideward direction, 
estimated values of Ri are greater than those of R s at each multiplicity n or h?. Both radii 
increase as multiplicity n increases. On the other hand, both radii decrease as kx increases. 

In the analysis of observed two-particle Bose-Einstein correlations at fixed multiplicities, 
the MD calculated from Eq.( l23|) becomes broader as multiplicity n increases, whereas it 
becomes narrower as kx increases. The behavior of calculated multiplicity distributions 
suggests that the data selection conditions in the observed two-particle Bose-Einstein corre- 
lations would have non-negligible effects on the shape of observed multiplicity distributions. 
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If the MD, one-particle and two-particle momentum densities are constructed from the 
same data samples, or at least from the data in the same pseudo-rapidity range, values of 
parameters would be estimated more consistently. 
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